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Section I 

INTRODUCTION 


Prior to the launching of a Space Shuttle Orbiter vehicle, a ferrying operation 
across the southern states ©ay be carried out to move the Orbiter from assembly 
site to launch point or from landing point to assembly site* The proposed 
route of this operation (Figure 1) has seven terminals and a like number of 
alternates for route segments averaging about 350 miles. The terminals are 
subject to certain landing and takeoff constraints and the route segments are 
subject to altitude and headvind constraints* The objective of this study is to 
find the distribution in time of flight duration between Edwards AFB and 
Kennedy Spaceflight Center (KSC) under these constraints, using meteorological 
data for the four midseasonal months of January, April, July, and October. 

To determine this distribution, two distinct statistical techniques have 
been used. The first is a Markov chain process and the second is a Monte Carlo 
method. A certain ferry configuration, a vehicle cruising speed of 250 knots, and 
some environmental constraints have been set for this particular problem. However, 
the two statistical methods are entirely applicable to other ferry configurations 
and other operational values of the variables pertaining to speed limitations and 
environmental requirements. 

To the best of the authors 9 knowledge, there is no previous work which is 
directly applicable to this topic. However , the Markov chain theory has been used 
by Brodie [1] in solving the similar problem of finding the probability of success 
of a space mission. A large amount of work has been done on model simulation with 
the Monte Carlo method (e.g., Schreider 12]). 

The study is divided into seven sections. The Markov chain theory is 
presented in Section II, with application to the simulation of Orbiter ferry 
flights. The flights are subject to certain persistence effects, so the Markov 
procedure is outlined for meteorological data which contain persistence as well 
as for data which are free of persistence. 
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Section III presents an adaptation of the Monte Carlo theory to the 
ptoble&t at hand, again treating the procedure for data containing persistence 
effects separately from persistence-free data* 

Section IV deals with the transition probabilities of Markov chain theory. 
After discussing the theoretical aspects of this topic, the application to ferry 
simulation is made by quantifying the Shuttle Orbiter ferry requirements* These 
constraints are related to hydrometeors, ceiling, visibility, inflight headwinds, 
and runway density. The effects of persistence in cruise headwinds and in 
hvdrometeors, or inclement weather, are also evaluated. 

Section V presents a number of distributions of flight duration under 
various assumptions, using both the Markov chain and Monte Carlo procedures. 

The degree of sensitivity of the results to changes in the cruise headwinds 
requirements and ceiling/visibility constraints is also tested, and the 
outcome when there is mutual dependence among the hydrometeors and the ceiling/ 
visibility constraint is investigated to give an optimistic, limiting assessment. 

Finally, the conclusions of this investigation are presented in Section VI, 
the references are shown in Section VII. 
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Section II 

MARKOV CHAIN THEORY 


Two applications of Markov Chain theory will be considered. The first 
case excludes the persistence of unfavorable conditions. The second case 
includes such effects. The two cases are presented in subsections 2.1 and 2.2, 
respectively. In subsection 2.3, the effect of additional independent delays is 
investigated. 


2.1 MARKOV THEORY WITHOUT PERSISTENCE EFFECTS 

Markov theory ueals with processes in which the probability of obtaining 

a particular state in the trial depends only upon the state preceding the 

trial. To illustrate, the probability of obtaining a chain of events E q , E^, ... 

E can he written as P<E , E, , ... E) * P(E ) P(E,/E ) ... P(E /E Here 

n oluolo n n~l 

P(E o > is the probability of being in the initial state E q . PtE^/E^) is the 
conditional probability of passing from state to in the first trial. 

In the simplest Shuttle ferrying case, the possible states can be assumed 

to be 


* Uj]. i* 1, 2* ...» 7 


( 1 ) 


where is the initial airport, a^ is the terminal for the first leg, etc. 
(see Figure 1), and &y is the final destination. Taking initially the case 
where all flights begin at the s<uae airport, then P(E o ) becomes 





1 


*<a 2 > 


0 

0 

♦ 

• 

* 

. 

• 


• 

P(a 7 ) 


» 

0 


m i * 1, 2, ... 7 


( 2 ) 


where [1^] is equivalent to P(E o >. 
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The term P(E /E , ) is the conditional probability of going from state E 
n n-l q- 

to state £ in the n th trial. This can be written as the matrix of transition 
n 

probabilities. 


P(E/E ) 
n ni 


1,1 

0 

0 

0 

0 

0 

0 


l U2 

r 2,2 

0 

0 

0 

0 

0 


2,3 

r 3,3 

0 

0 

0 

0 


0 

0 

r 3,4 

f 4,4 

0 

0 

0 


0 

0 

0 

r 4,5 

r 5,5 

0 


0 

0 

0 

0 

*5,6 

r 6,6 

0 


0 

0 

0 

0 

0 

r 6,7 

1 


(3) 


tT i,j ] ’ 1 


1, 2, • •*, j * 1» 2, • •*, 7 


Here, ^ is set equal to 1 because state is the absorbing state, or final 
destination* The zeroes represent the transition probabilities for nonadjacent 
airports; all flights operate between adjacent airports. 


If persistence effects are excluded, then the states E (or a^) refer to the 
airports along the route A^ (Figure 2a). Each trial is corridered to require half 
a day. The transition probability, gives the probability of leaving airport 

and reaching the next airport, A^ + ^, during the half-day trial period. The 
transition probabilities for morning trials are different from those for 
afternoon trials. 


In Figure 2a the initial state is shown at 1 (A.M.). This ?an be assumed 
to be N identical aircraft at the first airport ready for the first adorning trial* 

A fraction ^ of N will be successful and proceed as shown in Figure 2a to 
point a^, the second airport, at 1 (P*M.), ready for the afternoon trial* The 
remainder Nd-T^ * N ^ will be unsuccessful cad remain at a^ at 1 (P.M*) 
ready for the afternoon trial* 

After the afternoon trial, there is a probability T- * of the aircraft at 
airport moving to airport at 2 (A.K*)* 

The probability of reaching state a^, which is the final destination after 
a given number of trials, is the result being sought* The probability of attaining 
state ay after the first event is 
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P(l/a 7 ) = _ . . 

' 1-1 1 1 » / 


I.T, , • 0 


(4a) 


after 2 events. 


P(2/a,) 


7/7 \ 7 

:• I.j l T. T. ,\= • 1. T. . 

i=l \j'l i,7 J i*l 1 1,1 


( 2 ) 


(4b) 


The probability of reaching the final airport after n half-days is 


Z Z 7 

?(n/a 7 ) « • > ) 1 T, , T. , 

' 1=1 j=l m^l 1 1,J 


T = 7 I. 1. 

a, 7 i<L i 


(n) 


where 


1,7 


( 2 ) _ 


(3) = V T T T T 

U1 j=l k-1 Ui j ’ k k ’ 7 


7 7 


T. , T, - 
i.j j,7 


(°) m y •••y t t . . . t 
1,7 j»l m=l *->3 ^> k 31,7 


(4c) 

(5a) 

(5b) 

(5c) 


Since all aircraft are at the first airport initially, Eqs. (S) becomes 
7 7 

( 6 ) 


«»/.,) - £ <l)T lf} t 3>k - T„ >7 - T 1>7 

This matrix multiplication can be carried out directly. 


(n) 


s t “i p <°^ a 7 ) 

c n=0 7 


(7) 


is also of interest* It gives the probability of reaching the final airport 
th 

on or before the t trial* 


■r 


2.2 MARKOV THEORY WITH PERSISTENCE EFFECTS 

In this section » the term "persistence" signifies the effect of unfavorable 
conditions of the previous trial on the transition probabilities of the next 
trial. Such conditions signify that the previous trial was a "No go" outcome* 
Persistence is taken into account by redefining the Markov states as 
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E 


( 8 ) 



where a now refers to being located at the initial airport under previously 
favorable conditions; A^(F) a^ now refers to being located at the initial airport 
under previously unfavorable conditions: A^(U) a^ refers to being located at the 
next airport with previous conditions favorable; (F) is similar to but with 
previous conditions unfavorable, A^ (U) and so forth. Then as in Eq. (2), taking 
initially the cases where all flights begin at the same initial airport, P(E^) is 
redefined as 


P(E o ) 


~P(aj) ~ 


*P(F> 

P(a 2 ) 


P(U) 

P(a 3 ) 


0 

p<V 

S 

0 

_ P(a 13>_ 


• 

♦ 

0 


MV 


(9) 


where P(F) is the probability of the previous weather being favorable while 
P(U) is the probability of the previous weather being unfavorable. 


As before, the N identical aircraft are started at airport 1. However, 
this time a certain fraction P(F) have previously favorable weather and are in 
state while the other fraction P(U) * 1 - P(F) are considered to have 
previously unfavorable weather and are in state (Figure 2b). 

The transition probability T n * in Figure 2b gives the fraction of aircraft 

l* J 

at the first airport with favorable previous weather a^ at 1 (A.M.) that 
successfully reach the second airport, or state a^, at 1 (P.M.) after the first 
morning trial. T t 9 gives the other unsuccessful fraction which remain at the 
first airport but now with the previous weather unfavorable (e.g.» state a^) 
at 1 (P.M.) after the morning trial. 
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Similarly, the transition probability T., , in Figure 2b gives the fraction 
of aircraft at the first airport with unfavorable previous vea* ■ * state 

at 1 (A.M. ), that successfully reaches the second airport, or * tate a^ # . 1 (P.M.), 
after the morning trial. gives the unsuccessful fraction which remain at 

the first airporL in state at 1 (P.M.) ai :er the first trial. 


The probability of reaching the final destination in n trials is now given 


by 


13 

13 

13 


1 * * * f 

i*l 

i-i 

n^l 


13 


<n) 

M3 


(10a) 


The cumulative distribution is now given, as in Eq. (9), by 


n*0 


(10b) 


2.3 ADDITIONAL INDEPENDENT DELAYS 

During the flights, there will be occasions when an independent, additional 

delay beyond the previously calculated delays may occur* For instance a large 

scale weather formation may extend the overall flight time. The probability of 

an additional delay of d days of a flight can be written as p^. Then the 

probability P.(t) of reaching the destination in t days can be written as 
d 

p.(t) - p(t)u - i p.i + i P p(t-d) ai) 

d a-i d d-i d 

where P(t) is the probability of completing the trip in t days without the 
independent additional delay and P(t-d) is the probability of completing the 
trip in t-d days without such a delay. 

The cumulative distribution as in Eq. 7 will now become 

v. ■ !/(« <i2> 
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Section III 

MONTE CARLO THEORY 


In this section the Monte Carlo analysis is discussed* Subsection 3.1 is 
the case without persistence effects while subsection 3.2 extends the analysis 
to include the effects of persistence. 

3.1 MONTH CARLO THEORY WITHOUT PERSISTENCE EFFECTS 

Model sampling can be used to evaluate the probable number of days needed 

to reach the final destination. The simulated aircraft is started at the initial 

airport and its new location is calculated after the first event by the transition 

probability, T where T- , + T, 0 * 1. As before, the two subscripts indicate 

i*a i,2 

the airports of origin and of termination, respectively. A random number wh : ch 

is uniformly distributed between 0 and 1 is generated. If it is less than T. 

i 

the simulated aircraft is moved to the next airport. Otherwise, the aircraft 
is held at its present location. Its new position is calculated similarly for 
the next event. This process is continued until the simulated aircraft reaches 
the final airport aft< r t even s. It is then scored in the t-category. This 
process is then repeated for N samples. The probability of reaching the final 
destination after t trials can then be estimated by 

P(t/a ? ) * n t /N, (13) 

where n t is the number of simulations that reach the final destination in t 
trials, and N is the total number of simulations. 

The 95-percent confidence limits n n^/N, 6, are given by Schreider [2] as 

|n t /N - P(t/a 7 )| - i < 1.96 ^(t/a 7 )(l - P(t/a ? )/N - 1.96 o (14) 

This expression can be evaluated directly as before. 

3.2 MONTE CARLO THEORY WITH PERSISTENCE EFFECTS 

The Monte Carlo procedure can be modified, as was the Markov chain method, to 
include the effects of persistence. If P(F) is the probability of favorable 
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condition* «*" n given day* the simulated aircraft is . tried at the initial airport 

by picking a random number uniformly distributed i-ave;; t and l and tioparing it 

with P(F). If the random number is less than PfF), then the previous weather is 

considered favorable; otherwise, it is considered unfavorable* It the previous 

weather is favorable, i is used in the tiert trial; if not, T. is used. As 

l) i » 

discussed previously, this procedure is continued to find 

P(t/a l3 ) n t /N (1>) 




l 

I 



r i 

i I 
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Section IV 

TRANSITION PROBABILITIES 


In this section, the transition probal .ities used in the analysis are 
discussed. In subsection 4.1, the theoretical aspects are considered, first 
excluding persistence and then including this factor in the discussion. In 
subsection 4.2, the Shuttle Orbiter Ferry requirement * are quantified with 
persistence effects excluded. In subsection 4.3, the persistence effects are 
then added. 

4. 1 THEORETICAL ASPECTS 


* 


i 


4.1.1 Theoretical Aspects Without Persistence 

Consider a number of factors called A^, A^, ... A^. If any of these events 
occur, the aircraft will not proceed. The probability of events ot or ... 
or occurring is given by 


or A or ... or A ) = £ P(A ) - £ P(A.\) 

3=1 J j»k J k 


+ ^ ?(A 1 \ V ~ 

j,k,m 1 * m 

jfkfm 




j*k 

1 




(16) 


If A, 


and Aj are mutually exclusive. 


then 


P^A ) - P(A t ) P^/Aj) * 0 


(17a) 


* 



j ? 


* 




If A^ and A^ are iuctependent events, then 

P(AjAj) - P(A t ) P(A j /A i ) - P(A ± ) P(Aj) (17b) 

If A^ and are inclusive, so that whenever A^ occurs A^ must occur, 

P(A^Aj) * P(Aj/Aj) P(A t ) * ?(A ± ) (17c) 

These relationships will be used to evaluate the overall transition prob- 
abilities for the various operational constraints. 
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4,1*2 Theoretical Aspects With Persistence 

Consider as beiou- a nuab .r ;*r factors V, » A T . A r . tf anv or these 
events occur, the trial will be unfavorable, th**t 


» t * (A x or A, cr ...) £ 


(18) 


The probabt 1 ity of unfavorable wt-ai-ber in the present trial, given that the 
previous trial was unfavorable, is 

P(U t /U t-l ) = Pl[A J ° r A 2 or '*• Wl 1 

5 fk 


- I pua.) /u , i - ;; pf (a ,a > /u .] + 

1 i t t l j<k i k t t-l 

+ ji,. p “iV.> t , Vi> - - 


(19) 


* Tor T. , or ... 

2,2 4,4 

These quantities define the persistence factors desired in the calculation. 


If Aj is independent of the previous unfavorable weather, then 


Ft<AW Vl ] = PvA^) 


( 20 ) 


It can be seen that 


PtP t /0 t .l) • l - »tWl> - h.3 “ T 4,5 


( 21 ) 


which is the probability of a ”Go", or favorable weather on the present trial, 
given the previous trial was unfavorable, 0 


As given in Reference [3], 


P[U] P[F /U J 

Pl V F t-l 1 * 1 [1 - P(0)] * *1,3 ° r *3,5 ° r *** 


( 22 ) 
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From this result 


H[U t / Vl ] - I - HF t /F t .,] - I li2 .. t 3>4 or ... (23) 

4.2 QUANTIFICATION OF SHUTTLE ORBIT ER FERRY REQUIREMENTS 

The takeoff, landing, and inflight constraints specified for this study 
include: (1) operation under visual flight rules with no icing; (2) acceptable 

ceiling and visibility at the point of origin and the destination upon landing or 
taking off; (3) tolerable inflight headwinds; and (4) acceptable runway 
atmospheric density for takeoff. The first three of these constraints are 
somewhat interrelated, as is indicated in tfte Venn diagrams for eastbound and 
westbound flights (Figure 3). The aiount of interaction displayed in the 
diagrams is merely figurative. The symbols are explained in subsections 4.2.1 to 
4.2.4. 


To quantify the constraints, certain meteorological variables observed at 
the seven air bases were selected and their monthly summaries [4] were obtained 
over as long a period of record as possible. In general, the period of record 
exceeds 20 years. The meteorological variables which enter into the computation 
of transitional probability are listed in Table 4-1 with their specified con- 
straining values in some cases. The effects of possible long-term weather trends 
were not considered because of insufficient data. 

Table 4-1. SUMMARY OF THE METEOROLOGICAL VARIABLES UNOER FOUR CATEGORIES 
OF CONSTRAINTS, WITH AN INDICATION OF THE CONSTRAINING VALUES 


CATEGORY 

VARIABLE 

CONSTRAINT 

1 

Thunderstorm 

No operation 


Rain and/or drizzle 

No operation 


Freezing rain 

No operation 


Snow ••'n*7or sleet 

No operation 

2 

Cel lint/ visibility 

£1000 feet/3 miles 

3 

Headwinds 

Maxlirum wind velocity 

4 

Runway density 

Runway temperature < 103°F 


4.2.1 Constraints Rolotod to Hydromotoar* 

The variables under Category (1) are mutually exclusive in the sense that 
only one can be reported as "Present Weather" at a particular hour. Therefore, 
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the percentage frequency oi occurrence 01 the four variables in tnis category are 
simply added to form a hydrometeor group called ”H" in ti.«. Venn diagrams. These 
data are in the form of monthly summaries for each of eight 3 -hour periods per 
day. 

4.2.2 Constraints Related to Ceiling ond Visibility 

The variables under Category *.!} are reported jointly in the monthly 
summaries as the "percentage frequency of joint occurrence”. The probability 
of C/V for "No Go” due to this constraint was similarly taken from the monthly 
summaries. Thus the two most important parameters, ceiling and visibility, are 
conveniently combined in the available statistics. They are indicated as C/V 
in the Venn diagrams. 

4.2.3 Constraints Related to Inflight Headwinds 

The variable under Category (3) is taken into account by the use of 
equivalent headwind data supplied by the National Climatic Center, NOAA, 
Asheville, N. C. These data give the mean and standard deviation of the 
equivalent headwind for each of the six legs (Figure 1) and for each month. 

The 300-mb, 500-mb, and 700-mb levels are provided. 

To arrive at a transition probability for headwinds, a flight altitude 
of 3000 meters was selected for all legs except the westernmost leg between 
Edwards AFB and Tuscon. The latter segment was evaluated at 4600 meters. The 
specified constraints at these two levels are 9 meters per second and 14 meters 
per second, respectively. For east bound flights, a map inspection of several 
years 1 upper air charts showed that the probability of encountering headwinds 
greater than these magnitudes is negligible in all seasons. For westbound 
flights, a Gaussian distribution of equivalent headwind values for the given 
mean and variance is assumed, and the transition probability was calculated from 
standard tables. The headwinds parameter is called f V* in the Venn diagrams. 

4.2.4 Coaitrainti Rotated to Rimway Density 

The variable under Category (4) has probability values which are obtained 
from the percentage frequency of occurrence of temperature exceeding 100*F 
at the instrument shelter. This constraining value is believed to be a 
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conservative approximation to a corresponding "No Go" condition over the runway 
(currently assumed to be a runway temperature exceeding 103°F (see Table 4-1)). 
This variable is called "T" in the Venn diagrams. 

4.2.5 Transition Probability Volaos 

The transition probabilities corresponding to the four Categories are com- 
bined by the use of Eq. 16 (assuming terms beyond the second to be negligible) 
and Eq. 17b, yielding a set of values for each leg, for both eastbound and west- 
bound flights, for each of four midseasonal months (January, April, July, and 
October). These values, which are equivalent to are given in Table 4-2. 

For example, the value .0*7 which appears opposite 1000 hours under leg A is 
the probability of "No Go" in the morning from Edwards AFB to Tuscon in January. 

4.3 INCLUSION OF PERSISTENCE EFFECTS 

Tlie principal causes of delay due to persistence are found in (1) the 
winds at cruising levels and (2) quasi-s tat ionary or slowly moving systems 
breeding inclement operational weather over route segments and terminals. The 
first effect influences westbound flights because easterly winds are negligible; 
the second effect retards eastbound flights more than westbound flights because 
the motion of the synoptic systems are basically eastward and a flight usually 
cannot penetrate the system as it moves along its flight path. The essential 
information for counting consecutive days of delay, or "runs", was gleaned from 
an examination of eight years of Daily Weather Maps of the ESSA and NOAA organi- 
zations [5,6] subsequent to 1964. These Maps contain the surface and 500-mb 
charts, and prior to 1969 they include a second surface map spaced 12 hours from 
the max. chart. 

4.3.1 Inflight Hnodwiitd Persistence 

Headwind persistence was evaluated by noting the runs in which the opposing 
component of the wind at 500 mb was 40 kt or more* It is believed that this wind 
intensity at 500 mb usually indicates marginal conditions at lower cruising 
levels when the assigned constraints are observed. Three sectors consisting of 
two legs apiece were set up for the counting process. The result for each 
midseasonal month is given in Table 4-3. 
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Table 4-2. TRANSITION PROBABILITIES, Ti,i FOR EACH LEG OF EASTBOUNO 
AND WESTBOUND FLIGHTS IN EACH OF FOUR MIDSEASONAL MONTHS. 
THE EASTBOUND LEGS ARE DESIGNATED BY A (EDWARDS/TUSCON), 
B (TUSCON/EL PASO), C (EL PASO/ABILENE), D (ABILENE/ 
SHREVEPORT), E { SHREVEPORT / EGL IN), AND F (EGLIN/KSC). 

THE WESTBOUND LEGS ARE DESIGNATED BY A (KSC/EGLIN), ... 

F (TUSCON/EDWARDS). THE FLIGHT DEPARTURE TIMES ARE 
1000 AND 1300 LOCAL TIME 


EASTBOUND 

JANUARY 


A 

B 

C 

D 

E 

zrj 



T l,l 

■H 

EH 

ERE 

8ESPSI 

warn 

EH 

1000 

.047 


.193 

.334 

mmm 

.233 

1 

1300 

.043 


.151 

.493 

■s 

.185 

i 

1 

APRIL | 

1000 

.015 

.025 

.120 

■B 

.160 

.134 

1 

1300 

.030 

.041 

.098 

KB 

.140 

.134 

1 

JULY 

1000 

.002 

.018 

.053 

.067 

.160 

.160 

1 

1300 

.491 

.448 

.167 

.226 

.276 

.239 


OCTOBER 

1000 

.019 

.022 

.131 

.125 

.125 

.103 

1 

1300 

.036 

.035 

.089 

.084 

.017 

,099 

1 

WESTBOUND 

JANUARY 


A 

B 

C 

D 

E 

HCHR 

| 

wgm 

ran 

B) 1 

EH 

ERE 

IS 

EMI 

m 

1000 

.631 

.702 

.724 

.645 

.456 

.243 

n 

1300 

.610 

.676 

.702 

.629 

.435 

.240 

i 

APRIL 

1000 

■SB 

.542 



.607 

.530 

.358 

.157 

i 

1300 

H 

.528 

.546 

.519 

.367 

.169 

i 

JULY 

1000 

.220 

.238 

.150 

.297 

.064 

.014 

i 

1300 

.290 

.305 

.198 

i 

.242 

.205 

.421 

i 

OCTOBER 

— 

1000 

.319 

.332 

.363 

.363 

.202 

.069 

i 

1300 

.306 

.313 

.334 

.332 

.208 

.085 

1 
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Table 4-3. NUMBER OF RUNS OF OCCURRENCES OF FLIGHT DELAYS DUE TO 
PERSISTENT HEADWINDS, FOR WESTBOUND TRIPS IN JANUARY , 
APRIL, AND OCTOBER. THE SECTORS ARE KSC/SHR (KENNEDY 
SPACEFLIGHT CENTER/ SHREVEPORT), SHR/ELP (SHREVEPORT/ 

EL PASO), AND ELP/EDW (EL PASO/EDWARDS AFB). THE DATA 
ARE TAKEN FROM THE DAILY WEATHER KAP(s) SERIES 








RUNS 

(DAYS) 




N 

SECTOR 

2 

3 

4 

5 

6 

7 

8 

9 

10 

n 

12 


(Years) 













January 

7 1/2 

KSC/SHR 

4 

4 _ 

3 

3 

2 


2 



l 

2 



SHR/ELP 

8 

3 

2 

1 

1 


1 

3 

1 





ELP/EDW 

8 

4 

3 





1 




April 

8 

KSC/SHR 

4 

2 

1 

1 

1 



1 






SHR/ELP 

6 

3 

2 

2 

3 


2 



2 




ELP/EDW 

11 

JL 


1 

2 

2 






October 

9 

KSC/SHR 

2 













SHR/ELP 

2 

l 

1 











ELP/EDW 

3 

4 

i 










The presence of persistence in these data can be investigated by a method 
described by Brooks and Carruthers [7, pp. 309-313]. If there is no persistence, 
a theoretical distribution of runs of occurrences is computed by successive 

k 

evaluations of the quantity Np q, k « 1, 2, ... n, where N is the number of days 
in the sample; p is the independent probability of a constraining headwind, and 
is obtained as discussed in subsection 4.2.3; q * 1 - p; k is the number of days 
in the run. Thus the expected number of runs of at least 1 day, 2 days, ... n 
days are found. Taking the difference of adjacent values then gives the expected 
number of runs of exactly 1 day, 2 days, ... n-1 days. 

Applying this technique to the January headwind runs, which have a value of 
p * 0.555, the theoretical and empirical frequency distributions are computed and 
are presented in Table 4-4, just for the KSC/SHR sector. 
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Table 4-4. THEORETICAL AND EMPIRICAL FREQUENCY DISTRIBUTIONS FOR 
THE KSC/SHR {KENNEDY SPACEFLIGHT CENTER/SHREVEPORT) 
SECTOR IN JANUARY. N « 232 DAYS AND p - 0.555 


Calculated Value 
k or more days 
k days 

k 

Np k q 

A(Np k q) 

2 3 4 5 6 7 8 9 10 11 12 

57 32 18 10 5 3 2 1 1 

25 14 8 5 2 1 1 

Observed value 
k or more days 
k days 


21 17 13 10 7 5 5 3 3 3 2 

4 4 3 3 2 2 1 2 


The higher frequency in the last row of values, compared to the first 
row, when fc _> 6 indicates that persistence is present in the data. The next 
step is to assume that the probability of unfavorable conditions depends on the 
previous conditions, let this probability be designated as this is 

the probability of not reaching the next airport, given that headwind conditions 
prevented flight progress on the previous trial. 


Following Reference 7, the probability of 2 or more unfavorable days is 
given by P^. Then 3 or more unfavorable days is given by (P 2t ) P ^t^t-1^* and 
4 or more by (P^) Sunmin S> this gives 


S “ P 2t [1 + P(U t /D t-l ) + p2( Wi> + 


2t 


i - royop 


(24) 


which can be rewritten, where and S are obtained from the observed values 
(Table 4-4), as 


1 <vw 


NP, 


2t 


NS 


1 - 


21 


(21 + 17 + 13 + 


+ 2 ) 


0.764 


(25) 


Again, following Reference 7, a set of theoretical frequencies can be 
computed using * 0.764. The frequency of k or more unfavorable 

days is given by N(P^ t )P k 2 wbere ®£t *** is 

differences between the values on the top line of Table 4-5 give the expected 
number of runs of k unfavorable days. These values are compared to the observed 
values of runs of k unfavorable days as given in Table 4-4. 


17 




Table 4-5. THEORETICAL AND EMPIRICAL FREQUENCY DISTRIBUTIONS FOR 
THE SAME SECTOR AND MONTH AS IN TABLE 4-4, BUT WITH 
P(U t / U t . 1 ) - °. 7 64 


Calculated Value 


2 

3 

4 

5 

6 

7 

8 

9 

10 

11 12 

k or more days 


21 

16 

12 

9 

7 

5 

4 

3 

2 

2 1 

k days 

:.((21)P k - 2 (U t /U t _ 1 )} 

5 

4 

3 

2 

2 

1 

1 

1 

0 

1 

Observed value 












k days 


4 

4 

3 

3 

2 

0 

2 

0 

0 

1 


The encouraging outcome of this trial over the KSC/SHR sector led to the 
adoption of the above procedure to compute *(U t /U ) for each of the three 
sectors for .January, April, and October data. The results are listed in Table 4-6. 

Table 4-6. VALUES OF P(U t /Ut-i) COMPARED WITH VALUES OF P{U) FOR 
THREE WESTBOUND SECTORS IN JANUARY, APRIL, AND OCTOBER. 

THE SECTORS ARE AS DESIGNATED IN TABLE 4-3. 


MONTH 

SECTOR 


P(U) 

January 

KSC/SHR 

SHR/ELP 

ELP/EDW 

0.764 

0.710 

0.515 

0.555 

0.592 

0.316 

April 

KSC/SHR 

SHR/ELP 

ELP/EDW 

0.643 

0.733 

0.553 

0.435 

0.480 

0.243 

October 

KSC/SHR 

SHR/ELP 

ELP/EDW 

0.428 

0.428 

0.246 

0.270 

0.128 


When the effect of headwind persistence is taken into account by 
Eqs. (21-23) and the data in Table 4-6, a new transition matrix is produced. 

This matrix includes the probabilities related to the other constraints (see 
Figure 3) and the newly derived headwind probabilities, which exceed the old 
headwind probabilities in sectors where persistence is effective. The resulting 
time distribution of days required for a complete trip has been evaluated by the 
Markov two-state chain technique. The results are presented in Section V. 
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4.3.2 Hydrometeor Persistence 

The second effect of persistence has been analyze*! by scanning the 
surface charts during the period when two charts per day were published in the 
Daily Weather Map series. The most common typo of weather delay phenomenon is 
found to he a frontal wave in the Gulf of Mexico, with slow-moving cold fronts 
also accounting for a number of delay cases. Considering just east pound flights, 
which suffer more delays then westbound flights, the number of cases attributable 
to persistent synoptic- sc ale phenomena has been categorized according to length 
of delay for the entire trip. The data do not justify a sector by sector analysis, 
as was done for headwinds. The results, which are quite subjective, are 
presented in Table 4-7. 


Table 4-7. SUMARY OF THE EFFECTS OF PERSISTENT WEATHER DELAYS 
ON EASTBOUND FLIGHTS IN JANUARY, APRIL, AND OCTOBER. 
N IS THE NUMBER OF DAYS INCLUDED IN THE SURVEY 



N 

(Days) 

TOTAL DELAY 
(Days) 

NO. CASES 

NO. CASES/N 

January 

140 

1 1/2 

4 

.029 


2 

1 

.007 



2 1/2 

1 

.007 



3 

3 

.022 



3 1/2 

2 

.014 

April 

104 

1 1/2 

1 

.009 


2 

1 

.009 



2 1/2 

2 

.019 

October 

129 

1 1/2 

1 

.008 



2 

1 

.008 

i 

i 


2 1/2 

1 

.008 


The wrath of January is the only month studied which accrues an appreciable 
number of persistent weather interruptions of schedule. The probability of 
occurrence of an independent delay of t days was defined in subsection 2,3 (Eq. 
U). In Table 4-7, the values of p. are (see righthand column) p- * .029 + .007 * 

U ^ 

,036, p^ * .007 + .022 * .029, * .014. Since no transition probabilities are 

calculable for this particular kind of delay, the results are deferred to 
subsection 5.3. 
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Section V 
RESULTS 


Using the procedures which have been outlined in Sections il to IV, a 
number of computations have been carried out for the designated Shuttle Orb iter 
Ferry Route in compliance with certain en*. ironmental constraints* The results 
of these computations are presented in the following order: 

1. Cumulative time distributions of flight duration obtained by the 
Markov chain method, omitting persistence effects* 

2. Time distributions of flight duration obtained by the Monte Carlo 
process, omitting persistence effects. 

3. Same as 1, with persistence of inflight headwinds included. 

4. Tim*- distribution of flight duration in January based upon a survey 

of meteorological charts, including persistence effects of hydrometeors* 

5. Sensitivity studies of the variation of constraints for inflight head- 
winds and ceiling/visibility, using the Markov chain technique without 
persistence effects. 

6. Sensitivity study of an assumed interdependence among the hydrometeors, 
ceiling, and visibility, using the Markov chain technique without 
persistence effects. 

S.1 MARKOV CHAIN CUMULATIVE TIME DISTRIBUTIONS 

The Markov chain procedure applied under the designated constraints dis- 
cussed in Section IV* and confined to takeoff times of 1000 and 1300 local time 
with no persistence effects operating, yields the cumulative time distributions 
of flight duration between Edwards AFB and KSC shown in Figure 4. This outcome 
is also the probability of reaching the final terminal on or before a given day. 
For example, in Figure 4a the July curve indicates that 87 percent of the flights 
in this month should traverse the eastbound route from Edwards AFB to KSC 
in 4 days or less. Alternatively, there is an 87-percent probability that a 
particular eastbound trip will complete its flight in 4 days or less. The 
principal difference between eastbound and westbound results, as well as the 
difference between the seasons, is accounted for by the headwinds factor. The 
greater duration of July flights eastbound, compared to the transition seasons, 
is caused by the occurrence of precipitation forms and high summer temperatures 
(in the Southwest). 
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5.4 ADDITIONAL INDEPENDENT DLt V. 1 

The application of Kq. (1 » the data of subsection 4.J and the additional 

independent delays, or p, - vain. »f subset* tien i.'3, Yields a new time distri- 

d 

but ion *>f trip duration for east - id flights In January* The January probabilities 
«»re plotted (Figure 5a) U. sh • elongation oi the tail of tut* disit Ibutl^n 

*t tin* expense of shorter it* -.at ions of 3, 4, and 5 davs. The resulting 
*- ar-.uj.tt i vo time dfstrih* « % also plotted (Figure 4a) to disclose an increase 
of a* ft- day :• ti-* * * t r.-ent probability level of trip duration* Tnis 
itureas" ;tt t r : : utas *.e . * the effects of persistently inclement weather. 

Such delays are sia-ia&l in other midseasonai months (Table 4-7) and they are 
appreciably less on westbound flights than easthomid trips. Therefore, they are 
not represented additionally in graphs or tables. 

5.5 SENSITIVITY STUDIES 

A number of computer runs were made to observe the outcome when certain 
Shuttle Orbiter Ferry requirements were eliminated or modified. First, the 
inflight headwinds constraint was studied with the other constraints held at their 
nominal values. Second, the ceiling/visibil ity constraint was studied with the 
other constraints eliminated. 

5.5.1 Inflight Headwinds 

The results of a step by step relaxation of the designated headwind con- 
straint is related here, the limiting headwind values of 9 m sec * at >000 m and 
14 m sec * at 4o0G m being multiplied by factors of 1.25, 1.50, 1.75. and 2.00. 

When the four midseasonai months are combined, the resulting set of cumulative 
percentage frequency curves show a marked improvement as the limitation is 
relaxed progressively (Figure 8). For example, curve "a" corresponding to the 
designated nominal maximum headwind indicates that 95 percent of all westbound 
trips traverse the KSC/Edwards AFB route in 10 days or less. However, if the 
headwind constraint is relaxed to twice its designated value, the result (seen m 
curve M e") indicates that 95 percent of all westbound trips traverse the route 
in 6 days or less* The conclusion reached by this sensitivity study is that the 
longer durations of westbound trips in winter and transition seasons can be 
reduced significantly if greater headwind strengths can be toierated. 
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shows 


The alternative mode of presentation of Markov **!tain rv‘*ults <i*'igut«* ; 
just the number of days requited for 95 percent tiouiiJei l«»n **J ! lights. In January, 
12 1/2 days are nested tor l >> percent o! the flight*? to travel this route vest- 
l>ound under the nominal headwind constraint- Only 7 .lavs art* needed in the 
latter case if the headwind constraint i relaxed to twice uie designated value- 

5-5.2 Ceiling Visibility 

A second exper >metit tests The sens! t i vit v if the most fmportaa* weather** 
related variables- namely* the ceiling and visibility par.tru u-i when other 
constraints arc nil. The designated ;anding and takeoff minima are 1000 ft/3 mi 
/Category ”c ” ) * If these limits are relaxed to 500 ft/1 ml (Cat**gor> "d”), the 
time distribution can be computed as before for comparison (figure 10). Two 
other sets of limits are added for a sensitivity die- at tne 93-pet cent level. 

When Categories "c" and “d" are compared, the midseasonal months of Junu. ry, 
April, and October show a gain of about 1/2 day while duly lias almost no gain- A 
tightening of constraints to 2000 ft/5 Bii (Category "b") yields a loss of nearly 
one day in January and a loss of about 1/2 day in the other nrouths. A further 
tightening of constraints to 3000 ft/b miles (Category "a”) results in an 
additional extension of merely 1/2 day to one day in each month. This test 
therefore reveals little sensitivity in the ceiling and visibility* in comparison 
with the headwind constraint* 

5.5 INTERRELATIONSHIPS AMONG HYDROMETEORS, CEILING, AND VISIBILITY 

Although the above results assume that the events (e.g*, the meteorological 
record of constraining factors) are independent, but not mutually exclusive, this 
assumption is not strictly defensible. The degree of correlation between the 
various factors is not known. However, some interdependence can be expected 
between the hydrometeors and the ceiling/visibtlity parameter. A very optimistic 
assumption would be to state that such events are inclusive in the sense that the 
greatest probability of occurrence within Categories 1 and 2 (Tabie 4-1) represents 
that pair of Categories* When this substitution is made, the new set of transition 
probabilities gives a more optimistic result which may be regarded as a lower 
bound on the number of days required for a complete trip in each season* The 
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outcome for east bound and westbound trips, using the nominal headwind limits 
specified for the problem, is entered in Figure 7. Thus, in January the 
optimistic limit for eastbound flights is 4 1/2 days, compared to 5 days for 
the nominal case, whereas this limiting value for westbound flights is 11 1/2 
days, compared to 12 1/2 days for nominal ceiling and visibility constraints* 
Thus, the reduction of trip duration attributable to interdependence among 
hydrometeors, ceiling, and visibility is rather minimal* 
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Section Vi 
CONCLUSIONS 


The Markov chain and Monte Carlo analyse* are effective methods to 
determine probable flight times for Space Shuttle ferrying oy*. rat ions. Further- 
more, the techniques used in this study are applicable to other routes and other 
transport configurations* 

because of the assumption that only one leg can be completed each half-day, 
the present results indicate that the designated Shuttle Orbit or route and 
requirements correspond to an absolute minimum period of 3 davs. This period 
is needed to traverse the six-leg path beiv.^u Hdwards AFb and Kennedy Space- 
flight Center, flying in cither direction. 

East bound flights are found to be free of ground delays caused by the 
designated headwind constraint of 9 m sec * at 3000 m and 14 tn sec 1 at 4600 m* 
However, the other constraints (e.g. , ceiling, visibility, and hydrometeors) 
result in a 95-percent probability of completion within 5 days in any season* 
Conversely, there is a 5-percent risk ol exceeding 5 days in a flight started 
at random. Persistently unfavorable weather conditions extend the expected 
delay period an additional day in winter. 

Westbound flights are affected by the designated headwind constraint in 
winter and in transition seasons. Evaluation of the ground delays caused by 
this factor and the other constraints reveal an annual 95-percent completion 
level ranging from 5 days in July to 13 days in January. Inclusion of persistence 
in the headwinds has little effect except in winter, when the 13-day figure is 
raised to 17 days. 

Sensitivity tests upon the headwind constraint indicate that relaxation 
of this factor significantly reduces the 95-percent probability time. For 
example, increasing the limit to 1 1/2 times its designated value reduces the 
13-day figure to 9 days. On the other hand, the adjustmem of ceiling/visibility 
constraints has little effect upon flight duration. 
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Figure!. ROUTE OF THE SPACE SHUTTLE ORBITER FERRY VEHICLE. THE TERMINALS 
NUMBERED 1 TO 7 HAVE THE FOLLOWING NAMES: 1 - EDWARDS AFB, 2 - 
TUSCON (DAVIS-MONTHAN AFB), 3 - EL PASO (BIGGS AFB), 4 - ABILENE 
(DYESS AFB), 5 - SHREVEPORT (BARKSDALE AFB), 6 - EGLIN AFB, 

7 - KENNEDY SPACEFLIGHT CENTER 
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Figure 3. SCHEMATIC VENN DIAGRAM OF THE CONSTRAINING FACTORS FOR EASTBOUND AND 
WESTBOUND FLIGHTS* SHOWING A DEGREE OF INTERRELATIONSHIP. THE 
LETTER DESIGNATORS HAVE THE FOLLOWING MEANING: C/V - CEILING/ 
VISIBILITY, H - HYDROMETEORS, T - RUNWAY TEMPERATURE, W - HEADWINDS 
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Figure 4. CURVES OF CUMULATIVE PROBABILITY FOR NIWBER OF DAYS (t) REQUIRED TO FLY 
BETWEEN EDWARDS AFB AND KSC, EASTBOUND (4a) AND WESTBOUND (4b), FOR 
JANUARY, APRIL, JULY, AND OCTOBER. THE CONSTRAINTS LISTED IN TABLE 1 
ARE APPLIED, WITHOUT PERSISTENCE EFFECTS, TO OBTAIN THE CONTINUOUS 
CURVES. THE EFFECTS OF PERSISTENT, INCLEMENT WEATHER (SECTION V) 

ARE INCLUDED TO OBTAIN THE BROKEN CURVE IN 4a. THE MARKOV CHAIN 
PROCEDURE IS THE RELEVANT TECHNIQUE 
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Figure 6. PROBABILITY OF THE NUMBER OF DAYS (t) REQUIRED TO FLY WESTBOUND TRIPS 
FROM KSC TO EOWARDS AFB IN JANUARY, (6a), APRIL (6b), AND JULY (6c). 

. THE DOTS AND SHORT HORIZONTAL LINES ARE THE VALUES AND TWO-SIGMA 

CONFIDENCE LIMITS, RESPECTIVELY, FOR CALCULATIONS MADE BY THE MONTE 
CARLO METHOD. CIRCLES REPRESENT VALUES FOUND BY THE MARKOV CHAIN 
. PROCESS. 
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Figure 8. ANKlAl CURVES OF CUMULATIVE PROBABILITY FOR NUMBER OF DAYS (t) REQUIREO 
TO FLY BETWEEN EDWARDS AFB AND KSC, EASTBOUND AND WESTBOUND, FORMED BY 
AVERAGING THE DATA FOR JANUARY, APRIL, JULY, AND OCTOBER. CURVES a 
(WESTBOUND) AND X (EASTBOUND) ARE BASED UPON THE. DESIGNATED HEADWIND 
CONSTRAINTS OF 9 m sec-l AT 3.0 km AND 14 m sec" 1 AT 4.6 km, THE TWO 
FLIGHT LEVELS USED IN THE COMPUTATIONS. CURVES b, c, d, AND e (ALL 
WESTBOUND) ARE BASED UPON A RELAXATION OF THE HEADWIND CONSTRAINT , IN 
STEPS OF ONE-QUARTER, TO A MAXIMUM RELAXATION OF 2.00H. H REFERS TO 
EITHER OF THE DESIGNATED HEADWIND CONSTRAINTS. THE MARKOV CHAIN PROCESS 
IS THE RELEVANT METHOD 
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CATEGORY CIG (FT)/VSBY (HI) 


d 


c 

CONSTRAINT 

CATEGORY 

b 


a 

3 4 5 6 

t (Days) 


' 3000/6 

2000/5 

c 1000/3 

d 500/1 



Figure 10. 95-PERCENT PROBABILITY VALUES FOR THE NUMBER OF DAYS (t) REQUIRED TO 
aY EASTBOUND OR WESTBOUND TRIPS IN JANUARY, APRIL. JULY, AND OCTOBER, 
UNDER FOUR COMBINATIONS OF THE CEILING/VISIBILITY CONSTRAINT. THE 
OTHER CONSTRAINTS ARE NIL. THE MARKOV CHAIN PROCESS IS THE RELEVANT 
METHOD 
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